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ABSTRACT 



O ■ We present 5-15 /im imaging spectroscopy of the reflection nebulae vdB 17 

O : (NGC 1333), vdB59 (NGC 2068), vdBlOl, vdBlll, vdB 133 and vdB 135, 

^^ ■ obtained with the infrared camera and circular variable filter wheel on the 

Infrared Space Observatory (ISO). These nebulae are illuminated by stars with 
Teff = 3,600 K - 19,000 K, implying ultraviolet (UV; A < 400 nm) to total stellar 
CN . flux ratios of F(^x<monm) / Ftotai = 0.01 - 0.87. We detect the infrared emission 

features (lEFs) at 6.2, 7.7, 8.6, 11.3, and 12.7 /im, broad emission features at 6 - 
9 //m and 11 - 13 /xm, and 5-15 yUm continuum emission, from the interstellar 
0\ ■ medium in vdB 17, vdB 59, and vdB 133 (-F(A<4oonm) /-^totai = 0.22 - 0.87), and 

place upper limits on the emission from the interstellar medium in vdB 101, 
vdB 111, and vdB 135 (F(A<4oonm) /^tota/ = 0.01 - 0.20). 

Our goal is to test predictions of models attributing the lEFs to polycyclic 
c/3 . aromatic hydrocarbons (PAHs). Interstellar models predict PAHs change from 

singly ionized to neutral as the UV intensity. Go, decreases. The ratio of PAH 



emission at 6 - 10 /im to PAH emission at 10 - 14 /im is expected to be ten 
times higher in ionized than neutral PAHs. 

We observe no spectroscopic differences with varying Teff. We analyze the 
spectra of vdB 17 and vdB 59 as a function of distance from the star, to see how 
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the spectra depend on Gq within each source. The only quantitative difference 
we find is a broadening of the 7.7 /im lEF at Gq = 20 - 60 within vdB 17. We 
observe only a 40% change in the 6-10 /im to 10 - 14 /im flux ratio over Gq = 
20 to 6 X lOl 

Subject headings: reflection nebulae — dust, extinction — ISM: lines and bands 
— infrared: ISM: lines and bands — infrared: ISM: continuum 



1. Introduction 

The infrared emission features (lEFs) at 3.3, 6.2, 7.7, 8.6, 11.3 and 12.7 /im, and their 
associated continuum, have been the subject of intense study over the last decade. The 
lEF strengths correlate with the C/0 ratio in planetary nebulae, suggesting a carbon-based 
carrier (Roche & Aitken 1986; Cohen et al. 1986, 1989). Duley & Wilhams (1981) were the 
flrst to comment that some lEF wavelengths are characteristic of the bending and stretching 
modes of various C-H and C-C bonds in aromatic hydrocarbons. There has since been 
a steady accumulation of evidence, both laboratory and observational, in support of this 
assertion. The speciflc nature of the carrier (s), however, is yet to be determined. The lEFs 
have been attributed to polycyclic aromatic hydrocarbon (PAH) molecules (Leger & Puget 
1984; Allamandola, Tielens, & Barker 1985) and to more amorphous materials containing 
aromatic hydrocarbons (Sakata et al. 1984, 1987; Borghesi, Bussoletti, & Colangeli 1987; 
Blanco, Bussoletti, & Colangeli 1988; Duley 1988; Papoular et al. 1989). Allamandola, 
Tielens, & Barker (1989) and Bregman et al. (1989) attribute the broad spectral structures 
at 6 - 9 fira and 11-13 /im, which accompany the lEFs, to PAH clusters and amorphous 
carbon grains. 

The origin of the 1-25 /xm continuum emission, like that of the associated lEF 
emission, is still the subject of debate. The continuum is observed in regions of very low 
radiation density, such as reflection nebulae and infrared cirrus, where the grain equilibrium 
temperature is too low to explain the emission as being thermal in nature (Sellgren 1984; 
Draine & Anderson 1985; Guillois et al. 1996). It has been attributed to single ultraviolet 
(UV) photon excitation of very small grains (Sellgren, Werner, & Dinerstein 1983; Sellgren 
1984), to PAH molecular fluorescence (Leger & Puget 1984; Allamandola et al. 1985, 1989; 
Puget & Leger 1989), and to hydrogenated amorphous carbon grain luminescence (Duley & 
Williams 1988; Duley 1988). 

The lEF carriers are a ubiquitous component of the interstellar medium (ISM). Past 
studies have found lEFs in emission near Galactic H II regions, reflection nebulae, planetary 
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nebulae, proto-planetary nebulae, and the diffuse ISM of our own and other galaxies (see 
reviews by Aitken 1981; Allamandola et al. 1989; Puget & Leger 1989; Geballe 1997; and 
Tokunaga 1997). The 6.2 /im lEF has recently been detected in absorption in the diffuse 
ISM (Schutte et al. 1996, 1998). The 3.3 /xm lEF may have also been detected in absorption 
in molecular clouds (Sellgren, Smith, & Brooke 1994; Sellgren et al. 1995; Brooke, Sellgren, 
& Smith 1996; Brooke, Sellgren & Geballe 1999). 

lEF spectra have been recently classified into several groups (Geballe 1997; Tokunaga 
1997). Class A spectra at 5 - 15 /xm display narrow lEFs at 6.2, 7.7, 8.6, 11.3 and 12.7 /im, 
as well as associated continuum emission. The 7.7 /xm lEF is invariably the strongest of the 
Class A narrow lEFs (Geballe 1997). By contrast. Class B spectra lack the 8.6 /im lEF and 
have notably broader features at 8 and 11.5 /im (Kwok, Volk, & Hrivnak 1989). Class B 
spectra sometimes also show a strong broad emission feature at 8.8 /im (Buss et al. 1993). 
The overwhelming majority of lEF spectra currently measured fall under the category of 
Class A. Sources with Class A lEFs are typically UV-rich sources, such as the diffuse ISM, 
reflection nebulae, H II regions, young stars, planetary nebulae and galaxies hosting star 
formation. By contrast, the half-dozen or so Class B spectra found to date arise from F and 
G-type post-AGB stars forming planetary nebulae (i.e. proto-planetary nebulae; Geballe et 
al. 1992). 

Well-developed theoretical work exists for the absorption and emission properties of 
PAHs, inspired by the hypothesis that PAHs are responsible for the lEFs. Mixtures of 
PAHs and other dust components have been used to model the UV and visible extinction 
curve and the infrared emission observed in the ISM (Puget, Leger, & Boulanger 1985; 
Chlewicki & Laureijs 1988; Desert, Boulanger & Puget 1990; Jobhn, Leger, & Martin 1992; 
Siebenmorgen & Kriigel 1992; Schutte, Tielens, & Allamandola 1993; Dwek et al. 1997; Li 
& Greenberg 1997; Silva et al. 1998). All of these models include a size distribution of 
PAHs. Most of these size distributions include large PAHs, expected in the ISM but not 
yet measured in the laboratory, by extrapolation from laboratory data on smaller PAHs. 
Interstellar PAHs are predicted to be mainly positively charged and neutral in strong UV 
fields, and mainly neutral and negatively charged when the UV field is weak (Bakes & 
Tielens 1994, 1998; Salama et al. 1996; Dartois and d'Hendecourt 1997). Some interstellar 
models (particularly Desert et al. 1990) include both neutral and ionized PAHs. 

These models for the absorption and emission from a size distribution of neutral and 
ionized PAHs predict that the fraction of the total interstellar extinction curve due to PAH 
absorption is much larger at UV wavelengths than at visible wavelengths. This is vividly 
illustrated in Figures 2 and 3 of Desert et al. (1990), Figures 1, 2, and 3 of Joblin et al. 
(1992), Figures lb, 26 and 46 of Siebenmorgen & Kriigel (1992), Figure 4 of Dwek et al. 
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(1997), and Figure 3 of Silva et al. (1998). If PAHs contribute more to the total interstellar 
extinction at UV than visible wavelengths, then PAHs should absorb and re-emit a larger 
fraction of the total absorption and re-emission by dust in UV-rich sources than in UV-poor 
sources. 

Reflection nebulae are valuable tools for the study of the Class A interstellar lEFs. 
The localized heating of an interstellar cloud by a nearby, optically visible, star allows 
us to study the lEFs over varying excitation conditions by observing reflection nebulae 
illuminated by stars of different effective temperature, Tcs- This localized heating also 
results in brighter infrared emission than would be observed in a similar cloud heated only 
by the diffuse interstellar radiation field. 

Sellgren, Luan, & Werner (1990) analyzed IRAS observations of reflection nebulae 
illuminated by stars with Tes = 3,000 K - 33,000 K. They observed that the ratio of 
nebular emission within the IRAS 12 /xm broadband filter, Az//j,(12 yum), to the total 
bolometric far-infrared nebular emission, Jboi(FIR), was independent oiT^s, for T^g = 5,000 
K - 33,000 K. The total infrared emission of a reflection nebula is due to starlight absorbed 
and re-radiated by the mixture of molecules and grains responsible for the total interstellar 
extinction curve. Their observation that Az//,^(12 yum)//boi(FIR) is constant over a wide 
range of Tcs implies that the material emitting in the IRAS 12 /im band has an absorption 
curve with a wavelength dependence similar to that of the total interstellar absorption 
curve over both UV and visible wavelengths. 

The Sellgren et al. (1990) observations, finding Az/J,^(12 yum)/Jboi(FIR) to be 
independent of T^s, appear to be in confiict with theoretical predictions that the IRAS 
12 /im band is dominated by PAH emission and that PAHs contribute more to the total 
interstellar extinction curve at UV than visible wavelengths (Desert et al. 1990; Joblin et al. 
1992; Siebenmorgen & Kriigel 1992; Schutte et al. 1993; Dwek et al. 1997; Li & Greenberg 
1997; Silva et al. 1998). This is particularly puzzling as both Desert et al. (1990) and 
Siebenmorgen & Kriigel (1992) explicitly fit a mixture of PAHs and other dust components 
to the observed extinction curves and infrared emission of two refiection nebulae observed 
by Sellgren et al. (1990): NGC 2023 (vdB 52; Teg = 22,000 K) and NGC 7023 (vdB 139; 
Teff = 17,000 K). 

The Sellgren et al. (1990) observations could be reconciled with theoretical predictions 
if other ISM components besides PAHs contribute to the IRAS 12 /im broad-band emission 
in refiection nebulae illuminated by cool stars. Infrared spectroscopy can test whether 
the fractional contribution of PAHs to the IRAS 12 /im broad-band emission of refiection 
nebulae depends on Teg. Sellgren, Werner, & Allamandola (1996) have searched for spectral 
differences in the near-infrared emission of reflection nebulae with T^s = 3,600 K - 33,000 K. 
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They find that the equivalent width of the 3.3 /im lEF in reflection nebulae is independent 
of Teflf, over Tgfj = 11,000 K - 22,000 K, but they did not have the sensitivity to detect 
reflection nebulae illuminated by cooler stars. 

This paper presents a spectroscopic imaging study to search for and detail interstellar 
lEF emission toward a sample of reflection nebulae illuminated by stars with Tgfj = 3,600 
K - 19,000 K (Table |I]). We obtained our observations with the Infrared Space Observatory 
(ISO) camera, ISOCAM (Kessler et al. 1996; C. Cesarsky et al. 1996), using the circular 
variable filter (CVF). We have reported our initial result, the detection of lEFs toward 
the UV-poor nebula vdB 133 (illuminated by two stars with Tes = 6,800 K and 12,000 
K), in Uchida, Sellgren, & Werner (1998; hereafter Paper I). The goal of our study is to 
determine whether there are quantitative differences between the excitation of the different 
lEFs and the continuum emission among nebulae with widely varying Teg. We selected 
reflection nebulae with previous IRAS detections of extended 12 /im emission (Sellgren 
et al. 1990). The unprecedented sensitivity of the cryogenically cooled ISO satellite to 
low surface brightness mid-infrared emission has allowed us to search for lEFs toward the 
reflection nebulae illuminated by stars later than B-type, where the lEFs are too faint for 
detection with ground-based or airborne telescopes at ambient temperature. 



2. Observations 

2.1. Posit ion- Svi^itched Data 

We obtained ISOCAM multi-wavelength images of vdB 101, vdB 111, vdB 133 and 
vdB135 between 1996 March 19 and 1998 February 16. We used a set of 9 - 10 CVF 
wavelengths and 1-2 narrowband filters toward each nebula, in four sequences. Each 
sequence contains observations at three to four wavelengths, as listed in Table |^. We 
observed all four sequences for each nebula, except for vdB 135, where no Sequence 3 
observations were made. While the filter observations within a given sequence were 
performed consecutively in time, the scheduling of the sequences themselves were not 
necessarily constrained so. Some sequences thus share common wavelength observations to 
provide a check for consistency with time. All observations were made with a 6" pixel"^ 
scale giving a total field of view of 3.'2 with the 32x32 pixel array. The CVF spectral 
resolution was R = A/AA = 40 and the narrowband filters at 4.5 /im (LWl) and 15 /xm 
(LW9) had bandpasses of AA = 1 yum and 2 yum, respectively (ISOCAM Team 1994). 

The observations were position-switched (in a 2-cycle sky — >■ source —> sky —^ source 
sequence), with equal times on the source and sky to allow for the removal of the background 
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and dark current contributions. The source observations were centered on the illuminating 
star(s) of each nebula and the sky positions were toward nearby regions deemed free, 
based on inspection of the IRAS 12 fira skyview images, of any significant emission above 
the Zodiacal component and the local cirrus background. Table ^ lists the positions we 
observed. The basic integration time per exposure varied between 0.28 s and 10 s, depending 
on the wavelength and brightness of the central star. The total number of exposures was 
set to attain a total on-source integration time per filter of ~200 s for all observations. 



2.2. Full CVF Wavelength Scans 

Spectral images of vdB 17 (NGC 1333), vdB59 (NGC 2068), vdBlOl, vdB 133 and 
vdB 135 were also taken during the supplemental observing period of ISO. These involved 
full CVF wavelength scans at 5.14 - 9.44 /xm (CVFl) and 9.33 - 15.1 /im (CVF2), taken in 
both the increasing and decreasing wavelength directions, toward each of the source and sky 
positions. The spectral resolution was R = 40. The observations were spaced in wavelength 
to provide a factor of 1 - 2 oversampling in the vdB 17 and vdB 59 spectra, and 2-4 
oversampling in those of the other sources. The basic exposure time used was 2.1s, with 
8 exposures per CVF step (totalling 17.6 s) toward the bright sources vdB 17 and vdB 59, 
and 11 exposures per CVF step (totalling 23.1 s) toward vdB 101, vdB 133 and vdB 135. 
The pixel scale was 6" pixel"^ for the vdB 17, vdB 59 and vdB 133 measurements, like that 
of the earlier position-switched observations, and 12" pixel"^ for vdB 101 and vdB 135. In 
this second set of observations, the "source" fields were offset from the central star so as 
not to contain the star, thus allowing us to use longer integration times without saturating 
the detector. The vdB 133 field was centered toward the peak IFF emission in the region, 
as inferred from the earlier position-switched ISO observations, and those of vdB 17 and 
vdB 59 were centered in the directions containing the strongest extended infrared emission 
in our unpublished near-infrared ground-based images of the regions. The vdB 101 and 
vdB 135 fields were centered on the same positions observed by Sellgren et al. (1990) in 
their IRAS observations. The sky positions were chosen by examination of the 12 /im IRAS 
skyflux images. Table 3 lists the source and sky positions. 



3. Data Reduction 

The images were reduced, from the basic (CISP) data product, in the IDL-based 
ISOCAM data reduction environment known as "ICE", developed at the Institut 
d'Astrophysique Spatiale (IAS). By utilizing the basic ICE modules we were able to more 
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closely monitor and interact with the various stages of processing (i.e., the deghtching of 
cosmic ray hits from the data, detector transient response rectification and fiat-fielding) 
than with the standard (more automated) CAM lA processing package. Indeed, special 
efforts were made in the restoration of the source fluxes; the description of the method 
follows. 



3.1. The ISO Camera Transient Response 

The ISOCAM detector array experiences a response lag and a flux-dependent transient 
response (or residual memory effect). In order to recover the actual source fluxes, I{t), 
models of the detector behavior have been developed by several of the institutes involved 
with ISOCAM. We have adopted the transient model developed at the IAS (Abergel et al. 
1996), 

M{t) = r I{t) + (1 - r) / I{t') exp -— — - (1) 



oo 



which has proven from past experience to work best in rectifying low-level flux like that 
expected of the diffuse nebular emission. The measured flux, M, at time, t, is the sum of 
the instantaneous response of the detector, r I{t), and the integrated contribution of the 
exponentially decaying flux history of the detector (Abergel et al 1996). The detector is 
characterized by two parameters, r and a, the exact behavior of which is not well-known 
but which likely varies somewhat from pixel to pixel and perhaps with time, over the long 
term (weeks to months). The time constant for the exponential decay, r = a/M{t), is 
a function both of the measured flux and of the detector parameter, a. In the best of 
circumstances, the flux history of the detector prior to the start of one's own observations is 
known and can then be used to rectify the first part of the flux time series. Since, however, 
the data prior to one's own are generally not available, a number of "padding" frames, all 
with an estimated constant flux value, are instead prepended. 

If the detector response were perfect (Figure ^a), the flux history of a 2-cycle chop 
(sky — * source —^ sky —>■ source) sequence, for a single pixel or an average of pixels in the 
array, should appear as a series of square functions with the data points falling on two 
plateaus at a common minimum and two plateaus at a common maximum; the entire 
pattern then repeats itself (but likely with some general flux offset with respect to the first) 
for each of the subsequent position-switched cycles taken with a different filter. Figure |I|6 
displays an observed flux time sequence in raw form. The effects of the detector's finite 
response time and flux memory are seen to cause significant deviations from the ideal 



square pattern. The individual source and sky segments of the flux sequence, which should 
otherwise have constant levels, are either temporally increasing or decreasing (depending on 
the history of the detected signal) because of the non-ideal detector behavior. 

The position-switched data lend themselves well to the determination of the detector 
parameters. By using the ideal square-wave pattern as the target in a chi-square 
minimization procedure, with r, a, and the "padding" constant as free parameters, the best 
values of these parameters that rectify the observed time sequences can be found. Indeed, a 
set of parameters was determined for each of the four observation sequences toward vdB 133 
from averages of those pixels containing the nebular emission. The four parameter sets 
ranged between r = 0.72 - 0.73 and a = 954 - 1603 ADU gain"^, with padding values of 9.2 
- 13.8 ADU gain^^, where ADU gain~^ are analog-to-digital units divided by the detector 
gain. The specific parameter sets were used to rectify the four corresponding vdB 133 
observation sequences, and their average (r = 0.72 and a = 1200 ADU gain^^) was used 
to rectify the data sequences of the other sources. The parameter values determined by us 
provided some improvement over the results based on the values adopted by Abergel et al. 
(1996), r = 0.63 and a = 1200 ADU gain~^, for observations toward p Oph. Figure |l|c 
displays a rectified flux time sequence, after the observed flux time sequence (Fig. |l]6) has 
been rectified using the appropriate parameter values. 

In contrast, it is not possible to deduce the detector parameters from the full CVF 
scans with the method described above, since these scans involve sequential wavelength 
observations toward a given position and thus do not possess a predictable pattern, much 
less the regular "square wave" pattern of the position-switched sequences. All full CVF 
scan sequences of our supplemental observations were thus rectified using the average of 
the parameters determined from the original vdB 133 position-switched sequences: r = 0.72 
and a = 1200 ADU gain'^ 

The source and sky images were rectified separately before being used to produce the 
final source — sky maps. In each case, the dark current image (Boulanger 1996) was 
subtracted from the maps to establish the absolute count levels required for the rectification 
procedure. It was discovered, however, that the odd-even noise pattern (ISOCAM Team 
1994) remained in the images after a simple subtraction of the dark and that multiplication 
of the dark field (with factors ranging over 1.03 - 1.13) was needed to eliminate this 
problem. A single dark multiplication factor was used for all observations of a given source. 
An uncertainty of ~0.03 ADU gain~^ in determining the dark factor value translates to an 
uncertainty of ~3 MJy sr~^ in the individual source and sky map levels. However, since the 
source image and the sky image of a given nebula each use the same dark field (with the 
same dark multiplicative factor) , any errors in absolute flux levels of the component fields 
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do not propagate into the source — sky maps. The level offset in the source and sky images 
may lead to differential errors by the rectification process that appear in the source — sky 
maps, however, on the order of 5% of the absolute sky level (§4). 



4. Results 

Figures ^ ^ and ^ show our rectified, sky-subtracted spectra of reflection nebulae from 
position-switched observations and CVF scan observations. Each reflection nebula spectrum 
is a spatial average over a subsection of the image containing the localized lEF emission. 
The specifics of the spatial subregions used to produce the spectra are given in Table ^. 
For the position-switched data, where the image contained the illuminating star, the star 
and the immediately adjacent region containing a spurious arc-like reflection feature from 
the star (typically pixels X,Y = [6-20, 11-24] of the array) were excluded from the spatial 
integration. 

The position-switched observations were the first performed in this study. The 11 unique 
wavelength measurements were chosen to measure the strengths of the anticipated Class A 
interstellar lEFs (6.2, 7.6, 8.62, 11.22 /im), their surrounding continua (4.5, 5.8, 9.5, 10.5, 
and 15.0 /xm) and broad spectral structures such as the 6-9 fim bump and the 11 - 13 fim 
plateau (7.0, 8.4, and 12.0 fiia). 

The CVF scan observations, which provide a complete CVF spectrum over 5.14 - 15.1 
/xm, were performed second. In nearly all the CVF scans, there is a strong curvature at 
the start of the upward and downward scans because of the detector's flux memory and 
the satellite having slewed across the bright central star. With such large flux deviations, 
especially with the case of vdB 17 and vdB59, the response at the start of each scan could 
not entirely be corrected by the aforementioned rectification procedures (§3.1). To avoid 
these transient artifacts, we discarded data shortward of 8.0 /xm for the upward scan and 
data longward of 12.0 /im for the downward scan. Our final sky-subtracted spectra are 
thus constructed from the downward scan at 5.1 - 8.0 /im, from an average of the upward 
and downward scans at 8.0 - 12.0 yum, and from the upward scan at 12.0 - 15.1 /im. 
The difference between the up and down scans is plotted over a slightly wider range of 
wavelengths for each CVF-scanned spectrum, to estimate the uncertainty in our CVF scan 
spectra. 

The 9-10 /im segment of the CVF scan spectra also contain a prominent artifact 
resulting from the steeply declining transmission functions of the two CVF filter wheels 
near their wavelength overlap region: the long wavelength extreme of the CVFl wheel 
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(5.14 - 9.44 /iin) and the short wavelength extreme of the CVF2 wheel (9.33 - 16.52 yum). 
Thus the 9-10 /im spectral region, while shown in some spectra, has been excluded from 
detailed analysis. 

The position-switched spectrum of vdB 133 (Fig ^a) displays the distinct signature of 
Class A interstellar lEFs at 6.2, 7.6, 8.62, and 11.22 /im. The overwhelming majority of lEF 
spectra observed fall within this class (Geballe 1997; Tokunaga 1997). The full CVF-scan 
spectrum (Fig. ^b, and Fig. 1 of Paper I) taken during the supplemental observing period 
strongly confirms, both in location and relative intensities, the lEFs found in the more 
coarsely sampled positioned-switched spectrum. 

Figure |^ displays the CVF-scanned spectra toward the reflection nebulae vdB 17 and 
vdB59. The spectra of both nebulae clearly display the classic 6.2, 7.7, 8.6, 11.3, and 12.7 
/im lEFs that characterize Class A type interstellar lEF spectra. Both these nebulae are 
illuminated by hot, UV-rich stars {Tes = 12,000 K and 19,000 K for vdB 17 and vdB 59, 
respectively) and are typical of those that display strong interstellar lEFs. 

Figure ^ presents 5-15 //m spectra of vdB 101, vdB 111, and vdB 135 obtained with 
ISOCAM. Both the position-switched spectra and the fully-sampled CVF spectra of each 
source are shown, with the exception of vdB 111, toward which only position-switched 
observations were made. The fully sampled spectra of vdB 101 and vdB 135 are generally 
consistent with their position-switched counterparts obtained toward slightly different 
positions (Table ^. The source and background levels of vdB 101, vdB 111 and vdB 135 
measured by ISOCAM at 12 /im are consistent with those of IRAS. The reflection nebulae 
vdB 101, vdB 111, and vdB 135 are all illuminated by cool, UV-poor stars (Tefj = 5,000 K, 
7,300 K, and 3,600 K, respectively). These three nebulae have very low 12 /im IRAS Hnxes, 
ranging between 0.3 and 0.8 MJy sr~^. No lEFs, or continuum emission, are apparent in 
any of the source — sky spectra above the noise level of the observations, ~1 MJy sr~^. 

The task of detecting lEFs toward vdB 101, vdB 111, and vdB 135, even with ISOCAM 
operating under optimal conditions, was understood to be a difficult one. The predicted 
sensitivity of ISOCAM (~0.19 MJy sr~^ with 200 s integration time per filter in the 
position-switched observations) should have been sufficient, in principle. In practice, our 
actual ISOCAM sensitivity limit was ~1 MJy sr~^ for our observations, due to the faintness 
of the three nebulae compared to the relatively high background emission from diffuse 
Galactic emission along the same line of sight. The larger than anticipated uncertainty is 
attributed primarily to the limitations of the rectification procedure to completely correct 
for the detector effects and the uncertainties in the dark-current subtraction procedure 
required for the rectification process (discussed in §3). To achieve the expected sensitivity 
of 0.19 MJy sr~^ on top of a typical background level of 20 MJy sr~^ would require that 
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our flux rectification procedure be accurate to ~1%. Comparison of tlie upwards and 
downwards CVF scans indicated tliat in practice tlie rectification is only accurate to ~5%. 



5. Discussion 

In the following, we present a quantitative analysis of the lEFs as a function of the 
UV radiation field. We explore the effect of both the hardness of the UV field (Teg of the 
illuminating star) and the absolute intensity of the UV field (Go) • We expand our analysis 
over a wider range of Teg and Gq by comparing our observed sample to published spectra of 
vdBlOe (p Oph; Boulanger et al. 1996), and vdB 139 (NGC7023; D. Cesarsky et al. 1996). 
We also examine spectra at different distances from the illuminating star, for vdB 17 and 
vdB 59. This examines the effect of constant T^s but varying Go- 



5.1. Quantifying the UV field tov^^ard reflection nebulae 

The motivation for our ISO observations of reflection nebulae is the observed lack of 
dependence on T^s of Auli,{12 /xm)//boi(FIR) in reflection nebulae (Sellgren et al. 1990; 
see §1 and Table |I|). This is puzzling because current dust models predict that the IRAS 
12 /^m band is dominated by PAH emission and that PAHs contribute more to the total 
interstellar extinction curve at UV than visible wavelengths (Desert et al. 1990; Joblin 
et al. 1992; Siebenmorgen & Kriigel 1992; Schutte et al. 1993; Dwek et al. 1997; Li & 
Greenberg 1997; Silva et al. 1998). Because we are interested in analyzing ISO spectra as 
a function of UV strength, we need to quantify the UV strength in our sources. We have 
parameterized the UV fleld in four ways. 

Our flrst method of characterizing the UV fleld is simply by T^s for the illuminating star 
of each reflection nebula (Table |l|). Hotter stars emit a larger fraction of their total energy 
in the UV. This is complicated, however, by the difference in Tes for the two stars which 
illuminate vdB 133 [T^s = 6,800 K and 12,000 K; the cooler star is more luminous). 

Another measure of the UV fleld is the ratio of the stellar flux emitted at UV wavelengths 
to the total stellar flux. This simplifles the treatment of the UV fleld in vdB 133, since this 
ratio equals the sum of the UV fluxes of the two stars divided by the sum of the total flux 
of the two stars. The relationship between Teg and the ratio of UV to total stellar flux 
is, however, deflned by the longest wavelength included in the range of UV wavelengths. 
We have explored two ranges of UV wavelengths, both inspired by previous theoretical, 
observational, and laboratory work on PAHs and the lEFs. 
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The low ratio of visible to UV absorption, for the distribution of PAH sizes included in 
published models of interstellar absorption and emission, occurs because large molecules in 
general do not have electronic transitions beyond a cutoff wavelength, Acutofr, proportional 
to the molecule's length (Piatt 1956). Values of Acutofr = 300 - 500 nm have been estimated 
for neutral PAHs containing roughly 40 carbon atoms (Crawford, Tielens, & AUamandola 
1985; AUamandola et al. 1989; Leger et al. 1989; Desert et al. 1990; Joblin et al. 1992; 
Schutte et al. 1993). Ionized PAHs absorb at longer wavelengths than their neutral 
counterparts (Salama et al. 1996 and references therein). Approximate formulae for the 
relation between Acutofr and PAH size have been published for neutral PAHs (Desert et al. 
1990; Schutte et al. 1993) and for ionized PAHs (Desert et al. 1990; Salama et al. 1996). 

Thus, as our second method of parameterizing the UV field, we define a UV wavelength 
range of 91 - 400 nm. In this case, the fraction of total stellar flux emitted at UV 
wavelengths is -F(A<4oonm)/-^tota«- We chose Acutofr = 400 nm for this definition because it 
corresponds to the value of Acutofr, averaged over a distribution of PAH sizes and charge, 
at which current dust models predict the contribution of PAHs to the total interstellar 
extinction curve becomes negligible (Desert et al. 1990; Joblin et al. 1992; Siebenmorgen & 
Kriigel 1992; Dwek et al. 1997; Silva et al. 1998). Note 400 nm falls just longward of the 
observational dividing line between UV and visible wavelengths. 

Our third method of parameterizing the UV field is to define a UV wavelength range 
of 91 - 240 nm. Many previously published theoretical and observational papers on the 
ISM characterize the ambient radiation field by Go, the strength of the 91 - 240 nm UV 
radiation field. Examples of the prevalence of Go as a measure of the UV energy density 
include observational and theoretical studies of the 12 /xm IRAS emission in our own and 
other galaxies (Draine & Anderson 1985; Ryter, Puget, & Perault 1987; Boulanger et al. 
1988; Desert et al. 1990; Helou, Ryter, & Soifer 1991; Rowan- Robinson 1992), analysis of 
the observed strengths of the lEFs from recent ISO and IRTS spectra (Boulanger et al. 
1996, 1998a, 1998b; Onaka 1999), and models exploring the ionization state of PAHs and 
their effect on the heating and chemistry of different regions of the ISM (AUamandola et al. 
1989; Bakes & Tielens 1994, 1998; Spaans et al. 1994; Wolfire et al. 1995; Salama et al. 
1996; Dartois & d'Hendecourt 1997). For a UV range of 91 - 240 nm, the fraction of total 
stellar flux emitted at UV wavelengths is -F(A<240nm)/-^totai- 

Figure || displays F(^x<4oonm) / Ftotai and F(^x<240nra) / Ftotai as a function of T^s, determined 
using two different stellar models. We use Kurucz (1979) stellar atmospheric models for Tgs 
> 5,500 K, with log{g) = 4.0. We find, however, that the results are not sensitive to the 
value of log(5') adopted. Since the Kurucz models do not extend to cooler temperatures, we 
use a blackbody (Planck) function to model stars with Tes < 5,500 K. Figure |^ shows that a 
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blackbody slightly overestimates -F(A<4oonm) /-^totai and F(^x<2mnm) / Ftotai for the coolest stars 
for which we have Kurucz models. Figure ^ and Table show that F(^x<AQOnm) / Ftotai drops 
by a factor of 10, and -F(A<240nm) /-^totai falls by over a factor of 100, as Tgfj decreases from 
21,000K to 5,000 K. 

Our fourth method of parameterizing the UV field incident on a particular nebular 
position is by Go (sometimes called x or u), the strength of the 91 - 240 nm UV radiation 
field in units of the local interstellar radiation field. A value of Gq = 1 corresponds to a 
UV radiation field of 1.6 x lO'^ erg s'^ cm'^ at A = 91 - 240 nm (Habing 1968). Table 
H presents the distance to each source and the projected angular separation between the 
illuminating stars and the observed nebular region used to derive Gq. A single value of Gq 
suffices for vdB 106, vdB 133, and vdB 139, where the observed spectrum was obtained 
over a narrow range of distances from the central star, and for vdB 101, vdB 111, and vdB 
135, which we did not detect with ISO spectroscopy. For vdB 17 and vdB 59, we were able 
to obtain spectra at a range of values of stellar distance and thus Gq. Table ^ gives the 
value of Go at a typical nebular position for these two nebulae; we discuss how the spectra 
of vdB 17 and vdB 59 depend on nebular position (and thus Go) in a later section. 

We summarize the UV field for our observed sources and published comparison sources, 
as measured by Teg, -F(A<4oonm) /-Ftoiai, -F(A<240nm) /-Ftoiai, and Go, in Tables |l| and |. The 
observed values of Az//,^(12 /xm)/Jboi(FIR) (Sellgren et al. 1990; see §1 and Table |l|) are 
independent of all these measures of the UV strength. 



5.2. PAH ionization state in reflection nebulae 

Another reason to quantify the UV field in refiection nebulae is to search for spectroscopic 
differences between neutral and ionized PAHs. Interstellar PAHs are predicted to be mainly 
positively charged and neutral in refiection nebulae illuminated by B stars, and mainly 
neutral and negatively charged in the diffuse interstellar medium (Bakes & Tielens 1994, 
1998; Salama et al. 1996; Dartois and d'Hendecourt 1997). Predicted differences between 
neutral and ionized PAHs include changes both in the central wavelengths and in the 
relative strengths of PAH features (De Frees et al. 1993; Szczepanski & Vala 1993; Hudgins, 
Sandford, & AUamandola 1994; Pauzat, Talbi, & Ellinger 1995, 1997; Langhoff 1996; Cook 
& Saykally 1998; AUamandola, Hudgins, & Sandford 1999; Hudgins & AUamandola 1999b). 

The first ionization potential of PAHs lies in the range 4-8 eV, with a tendency for the 
ionization potential to decrease with increasing PAH size (Bakes & Tielens 1994; Salama et 
al. 1996). Hudgins & AUamandola (1999a) conclude from the wavelength spacing of the 6.2 
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and 7.7 /iin lEFs observed in the ISM that these lEFs arise from ionized PAHs with 50 - 
80 carbon atoms. A PAH with 50 - 80 carbon atoms has a first ionization potential of 5.8 - 
6.2 eV, corresponding to an ionization wavelength of 200 - 210 nm (Bakes & Tielens 1994). 
The value of Go (with Acutoff = 240 nm) is thus a rough indicator of the intensity of stellar 
radiation capable of ionizing PAHs at any given nebular position. 

Theory predicts that the mean ionization state of interstellar PAHs is a function of Gq 
(Bakes & Tielens 1994, 1998; Salama et al. 1996; Dartois and d'Hendecourt 1997). If the 
spectral features we observe in reflection nebulae are due to PAHs, then we would expect 
to detect changes both in the central wavelengths and in the relative strengths of lEFs 
between nebulae with different values of Gq. One goal of our observations is to test this 
prediction. 



5.3. Spatially averaged reflection nebulae spectra 

Localized enhancements of lEF emission, at 6.2, 7.7, 8.6, 11.3, and 12.7 /im, were 
detected in the source — sky spectra in three of the six nebulae observed: vdB17, vdB59 
and vdB 133. These three nebulae are illuminated by the hotter stars in our sample, with 
Teff = 6,800 K - 19,000 K. Our ~1 MJy sr~i upper limits on the lEF emission in vdB 101, 
vdBlll, and vdB 135 are consistent with their 12 /xm IRAS suTia.ce brightnesses of 0.3 - 
0.8 MJy sr^^ (Sellgren et al. 1990). These faint nebulae have cooler illuminating stars, with 
Teff = 3,600 K - 7,300 K. 

Figure ^ displays the superposed spectra of vdB 17, vdB 59, and vdB 133, those nebulae 
we observe having lEF emission in excess of the adjacent sky emission. These spectra 
(from Figs. ^ and §) are spatial averages over the emission region in each reflection nebula, 
defined by visual inspection of the images. Each spectrum is scaled by its average intensity 
over the 12 /^m IRAS hand (A ~ 7.5 — 15 /xm). We choose this scaling to illustrate the 
lEF peak and continuum levels (in MJy sr~^), normalized to an IRAS 12 /im intensity of 1 
MJy sr~^, for typical reflection nebulae. Such illustration is useful for planning for future 
mid-infrared spectroscopy with SIRTF and SOFIA of sources detected by the IRAS all-sky 
survey. 

The reflection nebulae spectra in Figure P, including their continua, broad emission 
bumps, and narrow lEFs, are roughly similar in shape despite the widely varying range of 
conditions quantified in Tables |l| and ^. The normalized peak fluxes of the different lEFs 
vary by less than a factor of two from source to source. More importantly, these variations 
show no systematic dependence on Teff. No evidence is seen among these three spectra for 
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disappearance of lEFs, appearance of new spectral features, shifts in lEF wavelengths, or 
dramatic changes in the strengths of lEFs relative to the broad emission bumps or the 
continuum emission. 



5.3.1. Quantifying the lEF Features 

The quantitative fluxes of the narrow lEFs, broad emission bumps, and continuum are 
quite sensitive to how the spectrum is divided into separate spectral components. Figure 
1^ displays the spectrum of vdB 17, to illustrate the range of possibihties for dividing the 
emission between the different spectral components. 

Figure ^a shows Method 1, in which we simultaneously fit Gaussians to the five narrow 
lEFs (at 6.2, 7.7, 8.6, 11.3 and 12.7 /im) and to two broad emission features at ~7 and 12 
/im. The latter two features represent the 6-9 /im "bump" (Cohen et al. 1986; Bregman 
et al. 1989) and the 11-13 /im "plateau" (Cohen, Tielens, & AUamandola 1985). The 
continuum is approximated by a spline fit through the data points at ~5.4, 10.3, and 14.8 
/im. The central wavelengths, widths, and heights of the Gaussians were all free parameters 
in the fit. Boulanger et al. (1998a) have suggested that the lEFs are best fit by Lorentzians. 
We used Gaussians rather than Lorentzians for the lEFs because they provide a better 
spectral fit to the lEFs observed in NGC 7023 aX R = 1000 (Sellgren et al. 1999). 

Figure |^6 shows Method 2, in which we simultaneously fit Gaussians to the five narrow 
lEFs (at 6.2, 7.7, 8.6, 11.3 and 12.7 /im). The continuum is modeled by a spline fit to all 
spectral data points outside that of the narrow lEFs; we assume the broad emission features 
at 6 - 9 /im and at 11-13 /im are not distinct from the continuum but rather are simply 
undulations in the continuum level in this case. Again the central wavelengths, widths, and 
heights of the Gaussians are free parameters in the fit. 

Figure |^c shows Method 3, in which we modeled the continuum by a linear fit to the 
spectral data points at ~5.3 and 15.0 /im, the wavelength extremes of the observations. No 
broad emission features are assumed in this case. The fiuxes of the five narrow lEFs were 
found by directly integrating over the continuum-subtracted spectra, the method most used 
in previous studies. The integration limits adopted for the 6.2, 7.7, 8.6, 11.3 and 12.7 /im 
lines were (6.0 - 6.6 /xm), (7.2 - 8.1 /xm), (8.1 - 9.9 /im), (10.5 - 12.0 /xm), (12.0 - 13.5 /im), 
respectively, as indicated in Figure |^c. The limits of the 6.2, 7.7, 8.6, and 11.3 /im lEFs are 
those used by Lu (1998). For the 12.7 /im lEF, which Lu does not include in his analysis, 
we integrate over a wavelength range similar in extent (AA = 1.5 /im) to that used for the 
11.3 /im line. 
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Figure illustrates the sensitivity of the estimated strengths of lEFs, broad features, and 
continuum to whether Method 1, 2, or 3 is used to fit the spectrum. The Gaussian used to 
fit the 6-9 /im broad bump in Method 1 is centered near 6.7 /im, and as a result the narrow 
7.7 /im lEF and especially the 8.6 fim lEF appear stronger. The higher adopted continuum 
level at 6 - 9 yum in Method 2 is centered closer to 7.5 /xm, which results in a somewhat 
weaker 7.7 /im lEF and a much weaker 8.6 /im lEF. The choice of a linear continuum 
(fitted to points at 5.3 and 15.0 /xm) in Method 3 consistently produce the highest 6-9 
/xm feature strengths of among the methods. The choice of broad integration limits for the 
8.6 /im feature (8.1 - 9.9 ^m) result in strengths comparable to or exceeding that of the 6.2 
/xm feature. Method 2 and Method 3 represent two extremes for the treatment of the 6-9 
yUm region. 

Previous studies of the relative strengths of different lEFs vary in their choices of 
continuum placement and integration limits in wavelength; the choice of continuum by 
Cohen et al. (1986, 1989) and Lu (1998) in their analyses is similar to Method 2, while that 
of Jourdain de Muizon, d'Hendecourt, & Geballe (1990) and Zavagno, Cox, & Baluteau 
(1992) is more like that of Methods 1 and 3. The sensitivity of our feature strengths 
to continuum placement, especially, demonstrates how much care must be exercised in 
intercomparing the results of different studies. 



5.3.2. Comparison of Nebulae 

Table |^ presents our quantitative results for the strengths of the lEFs, broad emission 
features, and continuum in vdB17, vdB59, vdB 133 (Paper I), vdB 139 (NGC7023; D. 
Cesarsky et al. 1996), and vdB 106 (p Oph; Boulanger et al. 1996) derived by three different 
methods for defining the lEFs and continuum. The methods by which the continuum 
and lEFs are defined, and thus the feature ratios quoted, vary between studies; the three 
methods employed here are representative of those typically used. The individual feature 
strengths are given in units of integrated line flux (ergs s~^ cm~^ sr^^) divided by the total 
integrated 11.3 /xm line flux (in ergs s~^ cm~^ sr~^). 

Table ^ shows that the relative strengths we derive for the lEFs, broad emission features, 
and continuum emission are very similar for all the sources we analyze, within a given 
method for decomposing the spectrum. The relative lEF strengths for any given source, 
however, are very sensitive to the method chosen. Figure |] illustrates this with the 8.6 to 
11.3 /im lEF ratio, 1(8.6 /im)//(11.3 /im), plotted against Teg. We originally chose this 
ratio because it is predicted to be much stronger for ionized PAHs than for neutral PAHs 
(De Frees et al. 1993; Szczepanski & Vala 1993; Hudgins et al. 1994; Pauzat et al. 1995, 



-17- 



1997; Langhoff 1996; Cook & Saykally 1998; Allamandola et al. 1999; see §5.4.1). The value 
of Tgfj (Table |^) is a measure of the hardness of the UV radiation field, which is independent 
of any adopted value of Acutoff, and which should be similar for all nebular positions 
included in the spatially averaged spectrum of each source. Figure ^ shows, however, that 
systematic differences in the mean values of 7(8.6 /im)//(11.3 /xm), measured by Methods 
1,2, and 3, are much larger than the source-to-source scatter of 7(8.6 /im)//(11.3 fim) for 
any particular method. 



5.3.3. Convolution of the Spectra with Broad-hand Filters 

Much effort has gone into characterizing the mid-infrared emission of the ISM using 
broadband filters, which provide limited spectral information but higher sensitivity than 
narrowband filters or spectrometers. The entire sky has been surveyed in the IRAS 12 /xm 
filter, and many ISOCAM studies have been conducted with the LW2 and LW3 filters, 
centered at 6.75 and 15 /^m respectively. Note that the ISOCAM LWIO filter and the 
ISOPHOT 11.5 /im filter are designed to be the same as the IRAS 12 /xm filter. Each 
of these broad filter bandpasses covers a mixture of lEFs, broad emission features, and 
continuum emission from interstellar dust and molecules. Table ^ therefore presents the 
expected IRAS 12 /im source fluxes and ISOCAM LW2 source fluxes found by convolving 
the source spectra with the 12 /im IRAS hand transmission profile and LW2 transmission 
profile (Fig. ^d), respectively. 

We have also investigated what fraction of the fiux within these broadband filters comes 
from each spectral component in the refiection nebulae we have observed. Table |^ gives 
the relative contributions of the lEFs, broad emission features, and continuum to the 
broadband fiuxes measured with broadband filters in common use. We derived these results 
by convolving the ISOCAM CVF spectra with the transmission profile of the LW2 (6.75 
/im) filter on ISOCAM and the IRAS 12 /im filter. Again, we find the estimated percentage 
contributions of lEFs, broad features, and continuum to the IRAS 12 /im filter and the 
ISOCAM LW2 filter are remarkably consistent from nebula to nebula, within a given 
Method (1, 2, or 3), but the absolute values for any nebula are quite sensitive to which 
Method (1, 2, or 3) is used to fit the spectrum. The choice of the continua via Method 3 
gives the maximal values of the narrow lEF/total ratios, while those of Method 2 represent 
the lower extreme. For the IRAS 12 /im filter, the narrow lEFs contribute typically ~20% 
(Method 1), ~40% (Method 2), and ~60% (Method 3) of the total fiux depending on the 
method used to define the continuum. All methods imply that 40 - 80% of the IRAS 12 /im 
fiux is due to continuum emission and/or broad emission features, rather than the narrow 
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5.4. Spatial variations within individual reflection nebulae 

The spectra for vdB 17 and vdB 59 plotted in Figures ^ and |^, and analyzed in Figures |^ 
and H and Table P, have been integrated over a wide range of projected distance from the 
central star. In this next section we explore how the spectral features within a single source 
image depend on position for vdB 17 and vdB 59. 

We have divided the spectrophotometric images of vdB 17 and 59 into radial bins, 
characterizing each bin by the incident UV starlight, Gq. We scale the values of Go in 
Table ^ by the inverse square of the projected angular separation of each pixel from the 
illuminating star. Figure ^ shows the average spectra of different radial bins within vdB 17 
and vdB 59. The values of Gq in these bins are equally spaced in log((jo); they vary over 
Go = 20 - 4 X 10^ for vdB 17, and over Go = 200 - 6 x 10^ for vdB 59. 

Near-infrared long-slit spectroscopy of vdB 17 and 59 (Martini, Sellgren, & DePoy 1999), 
as well as our unpublished ground-based near-infrared imaging and the ISO images we 
analyze here, suggest that vdB 17 has a much simpler geometry than vdB 59. Thus our use 
of the projected distance between a pixel and the central star to estimate Go for each pixel 
is likely to be a better approximation for vdB 17 than for vdB 59. 

Figure ^ shows interesting variations in the spectra of vdB 17 as a function of nebular 
position (or Go). At first glance, it appears that the 8.6 /im lEF steadily becomes less 
prominent in the spectra at low values of Gq. A more careful examination, however, shows 
that the peak intensity of the 8.6 fira lEF, relative to the peak intensity of the 7.7 /im lEF, 
is relatively constant with Go- Instead, it seems that the apparent weakening of the 8.6 fim 
lEF at low Go is rather due to a broadening of the overlapping 7.7 /xm lEF. 

Figure |10| quantifies the change in the full- width at half-maximum (FWHM) of the 7.7 
/xm lEF with Go in vdB 17. We adopted the method of Boulanger et al. (1998a) to measure 
the FWHMs of different lEFs, by fitting a linear baseline to continuum points near 5 and 15 
/xm, and then fitting a blend of three Lorentzians to the 6.2, 7.7, and 8.6 fim lEFs. Figure 
p!0| clearly shows that the 7.7 fim lEF is systematically broader at low values of Go in vdB 
17. Figure ^ also suggests that perhaps the 6.2 /im lEF may also broaden with lower Go, 
but the evidence is not as compelling as it is for the 7.7 /xm lEF. No change in the FWHM 
of the 8.6 /xm lEF is observed within the uncertainties. 

The spectra of vdB 59, in contrast to those of vdB 17, do not show any strong systematic 



-19- 



dependence on Gq. The width of the 7.7 /iin lEF and the contrast of the 8.6 /im lEF 
against the red wing of the 7.7 /im lEF show httle change in Figure ^. This may be due to 
the more complex geometry of vdB 59, where emission from regions close to and far from 
the star appear to contribute to the same line of sight (Martini et al. 1999). Alternatively, 
this may be because Go in vdB 59 does not extend to the low values (Go = 20 - 60) where 
the most marked spectral changes are observed in vdB 17. 

The origin of the widths of the lEFs is quite controversial and this complicates any 
explanation of why the 7.7 /xm lEF width changes spatially within vdB 17. An increase in 
PAH temperature should broaden the lEFs, but that would also result in a shift in the 
lEF central wavelengths (Joblin et al. 1995). We, however, do not observe any such shift 
in the 7.7 /xm lEF central wavelength. An increase in particle size might broaden the lEFs 
(Boulanger et al. 1998). A change in the mix of the lEF emitters, all of which have slightly 
different wavelengths, could also lead to a change in the lEF widths (Cook & Saykally 1998; 
Le Coupanec et al. 1998; Papoular 1999). Other broadening mechanisms are discussed by 
the above authors; whichever mechanism accounts for the broadening of the 7.7 /im lEF 
that we observe in vdB 17 needs to be consistent with the much lower, or non-existent, 
broadening we observe in the 6.2 and 8.6 /xm lEFs, and with the constancy in the central 
wavelengths of all the observed lEFs. 

The peak strengths of the 6.2 and 8.6 /im lEFs, relative to the peak strength of the 7.7 
//m lEF, show little variation among spectra at different distances within vdB 17 and vdB 
59, over Go = 20 - 6 x 10^ (Fig. |). The peak intensity of the 8.6 fim lEF of vdB 17, 
relative to the peak intensity of the 7.7 fim lEF, has values of 0.53 - 0.56 in Figure |^. A 
similar constancy in peak intensity (relative to the peak 7.7 /xm lEF intensity) is observed 
for the 8.6 /xm lEF within vdB 59 and for the 6.2 /im lEF within both vdB 17 and vdB 59. 

Figure ^ also shows marginal evidence for a weak emission feature near 6.9 /im appearing 
in the spectra of one or two of the most distant radial bins within vdB 17. A weak 
interstellar emission feature at 6.9 yum was first discovered by Bregman et al. (1983) and 
further analyzed by Cohen et al. (1986). This emission feature can be confused with [Ar II] 
emission at 7.0 /im, and so can only be reliably detected at low spectral resolution in sources 
free of ionized hydrogen, such as reflection nebulae. The 6.9 /xm feature is not detected in 
spectra of radial bins closer to the star in vdB 17 (Fig. ^. The 6.9 /xm feature is also not 
detected in the spectra of any radial bin for vdB 59 (Fig. ^, or in the spatially integrated 
ISO spectra of vdB 17 or vdB 59 (Fig. |), vdB 133 (Fig. g Paper I), vdB 139 (NGC 7023; 
D. Cesarsky et al. 1996; Laureijs et al. 1996) or vdB 106 (p Oph; Boulanger et al. 1996). 
The 0-0 S(5) H2 line at 6.9 /xm is unlikely to be a contributor, as other H2 lines with similar 
predicted strengths, at 8.0, 9.7, and 12.3 pm, are not detected, and because the 1-0 S(l) H2 
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emission in our unpublished narrowband images of vdB 17 is concentrated near the star 
at higher Gq values. We are therefore cautious in claiming that we have detected 6.9 /xm 
emission in the outermost regions of vdB 17. 



5. 4-1- The 8.6 ^m emission feature 

Laboratory data and theoretical calculations for PAHs show that the 8.6 to 11.3 /im 
lEF flux ratio, 1(8.6 /im)//(11.3 /im), is much higher in ionized than neutral PAHs (De 
Frees et al. 1993; Szczepanski & Vala 1993; Hudgins et al. 1994; Pauzat et al. 1995, 1997; 
Langhoff 1996; Cook & Saykally 1998; AUamandola et al. 1999). The 8.6 /im and 11.3 yum 
bands in PAHs correspond to the in-plane bend and out-of-plane bend, respectively, of the 
same aromatic C-H bond. Thus 7(8.6 /im)//(11.3 /im) should be a useful measure of PAH 
ionization because it should be less sensitive to PAH size or PAH hydrogenation than the 
ratio of an aromatic C-C vibration to an aromatic C-H vibration. 

The mean PAH ionization is predicted to increase with Go (Bakes & Tielens 1994, 1998; 
Salama et al. 1996; Dartois and d'Hendecourt 1997). One test of the PAH hypothesis, 
therefore, is to measure 7(8.6 /im)//(11.3 /im) as a function of Gq. Joblin et al. (1996) 
have observed 7(8.6 /im)/7(11.3 /xm) to increase closer to the star in the reflection nebula 
NGC 1333/SVS-3. They interpret this as due to PAHs becoming increasingly ionized near 
the B-type illuminating star. 

Figures |^, ^ and Table ^, however, demonstrate how strongly the measured value of 
7(8.6 /im)/7(11.3 /im) in our data depends on the method used to separate the 8.6 /im lEF 
from other spectral components such as the 7.7 /im lEF and the 6-9 /im broad emission 
bump. A local baseline (Method 2) to measure the strength of the 8.6 /im lEF would lead 
us to conclude that the 8.6 /im lEF virtually disappears at low values of Gq in vdB 17, as 
can be seen qualitatively in Figure ^ The observed broadening of the 7.7 /xm lEF at low Gq 
in vdB 17 (Figure 0) will, however, cause an increase in strength in the wing of the 7.7 /xm 
lEF. Because the red wing of the 7.7 /im lEF deflnes the local continuum for the 8.6 /im 
lEF in Method 2, any broadening of the 7.7 /im lEF (at constant 7.7 /im peak intensity) 
naturally results in an apparent weakening in the 8.6 /im lEF strength measured by this 
method. When we treat the 7.7 and 8.6 /im lEFs as a blend of two symmetrical features 
lying above a much lower continuum level (Method 1), then we flnd that the 8.6 /im lEF 
shows little variation with Go in vdB 17. 

We conclude that 7(8.6 /xm)/7(11.3 /im) is extremely dependent on how the 8.6 /im lEF 
is measured (Figs. and ||; Table ^. Joblin et al. (1996) use a local basehne (Method 2) 
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in measuring the strength of the 8.6 /im lEF, because their ground-based spectra do not 
extend to short enough wavelengths to include the entire 7.7 /im lEF. The increase in the 
7.7 /im lEF width with decreasing Gq we observe in vdB 17 (Figs. |^ and 0) leads us to 
urge caution in using 7(8.6 /im)//(11.3 /im) to quantify PAH ionization. 



5.4-2. The ratio of 6 - 10 fim emission to 10 - I4. fifn emission 

The difficulty of quantifying the observed 8.6 /im lEF strength inspires us to search for 
a more robust way to test predictions of the PAH ionization state as a function of Gq. 
Laboratory data and theoretical calculations for PAHs show that the sum of all the PAH 
bands at 6 - 10 /im, relative to the sum of all the PAH bands at 10 - 14 /im, is at least an 
order of magnitude higher in ionized PAHs than in neutral PAHs (De Frees et al. 1993; 
Szczepanski & Vala 1993; Hudgins et al. 1994; Pauzat et al. 1995, 1997; Langhoff 1996; 
Cook & Saykally 1998; Allamandola et al. 1999). We therefore have measured 7(5.50-9.75 
/im)/J(10. 25-14.0 /im), the ratio of the integrated intensity at 5.50 - 9.75 /im to the 
integrated intensity at 10.25 - 14.0 /xm, in ISO spectra of reflection nebulae, as an alternate 
way to search for the spectroscopic signature of changing PAH ionization state. 

Figure jl^ presents our measurements of /(5. 50-9. 75 /im)//(10. 25-14.0 /xm) as a function 
of Go ill reflection nebulae. We subtracted a continuum from each spectrum before 
measuring 1(5.50-9.75 /im)//(10. 25-14.0 /im). The 5.50 - 9.75 /im continuum was defined 
by a linear fit to data points near 5.3 and 10.0 /im. The 10.25 - 14.0 /im continuum was 
defined by a linear fit to data points near 10.0 and 14.5 /im. We plot a single value of 
/(5.50-9.75 /im)//(10.25-14.0 /xm) for vdB 106 (p Oph; Boulanger et al. 1996), vdB 133 
(Fig. gb and Paper I), and vdB 139 (NGC 7023; D. Cesarsky et al. 1996). For vdB 17 and 
vdB 59, we plot /(5. 50-9. 75 /im)//(10. 25-14.0 /xm) at six values of Go, derived from the 
spectra of the radial bins shown in Figure |^. 



The remarkable conclusion from Figure |TT] is that /(5. 50-9. 75 /im)//(10. 25-14.0 /xm) 



changes by no more than 40% among sources with Go = 20 - 6 x 10^. This observational 
ratio should not depend at all on how the observed interstellar spectra are divided into 
separate narrow lEFs and broader emission components. We emphasize that /(5. 50-9. 75 
/xm)//(10. 25-14.0 /im) should be a more robust and sensitive measure of the PAH ionization 
state than J(8.6 /xm)//(11.3 /xm). 

Our observation that the mid-infrared spectra of vdB 17 and vdB 59 are very similar 
for Go = 20 - 6 x 10^ (Figs. ^ and |llD are in strong agreement with other recent 



work. Boulanger et al. (1998b) concludes that the lEF fiux ratios in ISO spectra of the 
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Chamaeleon cloud, the p Oph cloud, NGC 7023, and Ml 7 show no systematic dependence 
on Go) over Gq = 1 - 10^. Onaka (1999) finds that the lEF flux ratios in ISO spectra of the 
Carina nebula are independent of Gq, over Gq = 100 - 10^. 



6. Conclusions 

We have obtained 5-15 /im ISOCAM CVF spectroscopy of six reflection nebulae, in 
order to test predictions of interstellar models for the infrared emission features. Our 
sample includes a range of UV fleld hardness {T^^ = 3,600 K - 19,000 K for the illuminating 

stars; F(^x<momn) / Ftotai = 0.01 - 0.87) and UV fleld intensity (Gq = 20 - 6 x 10^). 

Our flrst goal was to search for spectroscopic differences as a function of Tcs for the 
illuminating stars. Detection of such spectroscopic differences could potentially reconcile 
12 //m broad-band /i^yl^ observations showing that Az//i^(12 /im)//boi(FIR) is independent 
of Teff, for reflection nebulae illuminated by stars T^s = 5,000 - 33,000 K (Sellgren et al. 
1990), with model predictions that the IRAS 12 /im emission is dominated by PAHs which 
contribute much more to the total interstellar extinction curve at UV wavelengths than at 
visible wavelengths (Desert et al. 1990; Joblin et al. 1992; Siebenmorgen & Kriigel 1992; 
Dwek et al. 1997; Silva et al. 1998). 

A second goal was to search for spectroscopic differences as a function of Gq. Interstellar 
models predict that more PAHs are positively charged when Gq increases (Bakes & 
Tielens 1994, 1998; Salama et al. 1996; Dartois and d'Hendecourt 1997). Laboratory data 
and theoretical calculations for PAHs show that 1(8.6 /im)//(11.3 fim) and 1(5.50-9.75 
/im)//(10. 25-14.0 fim) are both much higher in ionized than neutral PAHs (De Frees et al. 
1993; Szczepanski & Vala 1993; Hudgins et al. 1994; Pauzat et al. 1995, 1997; Langhoff 
1996; Cook & Saykally 1998; Allamandola et al. 1999). 

We detect interstellar lEFs at 6.2, 7.7, 8.6, 11.3, and 12.7 ;um, broad emission features 
at 6 - 9 fim. and 11-13 /im, and continuum emission at 5 - 15 fira in the reflection 
nebulae vdB 17, vdB 59, and vdB 133. These nebulae were the brightest in our sample and 
illuminated by hotter stars (-F'(A<4oonm) /-^totai = 0.22 - 0.87). The normalized peak fluxes 
of the different lEFs vary by less than a factor of two among these three nebulae, with no 
systematic dependence on Tefj . No evidence is seen among these spectra for disappearance 
of lEFs, appearance of new spectral features, shifts in lEF wavelengths, or dramatic changes 
in the strengths of lEFs relative to the broad emission bumps or the continuum emission. 

We do not detect any feature emission or continuum emission in vdB 101, vdB 111, and 
vdB 135. These nebulae were the faintest in our sample and were illuminated by cooler 
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stars (-F'(A<400nm)/-^tota« = 0.01 - 0.20). Their upper limits are consistent with their low 
measured 12 /im IRAS surface brightnesses. 

We compare spatially averaged spectra of vdB 17, vdB 59, and vdB 133 to published 
spectra of vdB 106 (p Oph; Boulanger et al. 1996) and vdB 139 (NGC 7023; D. Cesarsky 
et al. 1996). Qualitative comparison of spectra show all these spectra are very similar over 
a range of UV hardness {T^s = 6,800 - 22,000; F(A<4oonm) /i^totai = 0.22 - 0.92) and UV 
intensity {Gq = 40 - 1800). 

We find that quantifying the lEF fluxes is quite sensitive to the method by which the 
spectra are separated into lEFs, broad emission bumps, and continuum. For instance, 1(8.6 
/im)//(11.3 fim.) shows no dependence on T^s for any given method of defining the lEF 
fluxes, but 1(8.6 /im)//(11.3 ^m) shows strong systematic differences between the values 
derived by different methods. 

We also quantify the contributions of lEFs, broad emission bumps, and continuum to 
broad-band filters such as the IRAS 12 /im filter and the ISOCAM LW2 (6.75 fim) filter. 
Again, for any given method of distinguishing between these separate spectral components, 
the contributions of each component to the broad-band filters are independent of Tes, but 
the absolute value of the contribution depends strongly on the method for how the spectrum 
is divided among lEFs, broad emission bumps, and continuum. 

We divide the spectra of vdB 17 and vdB 59 into radial bins, with the UV field of each 
bin characterized by Go- We find that the peak intensities of the 6.2 and 8.6 /im lEFs, 
compared to the peak intensity of the 7.7 fira lEF, show no variation over Gq = 20 - 6 x 
10^ in vdB 17 and vdB 59. 

In the most distant regions of vdB 17, we observe that the 8.6 /im lEF appears weaker 
due to increased blending with the 7.7 /im lEF. We quantify the FWHM of the 7.7 /im 
feature, and find that it grows systematically broader at low Go in vdB 17. We conclude 
that the low contrast of the 8.6 /im lEF at low Gq in vdB 17 is the result of the increased 
width of the 7.7 /im lEF. This observed broadening of the 7.7 /im lEF, combined with the 
sensitivity of /(8.6 /im)//(11.3 fim.) to baseline placement for the 8.6 /xm lEF, suggests that 
caution should be used in inferring PAH ionization from 1(8.6 /im)//(11.3 /im). 

We suggest using /(5. 50-9. 75 //m)//(10. 25-14.0 /im) as an alternate method of testing 
the predictions of PAH ionization as a function of Gq. We measure this in single nebular 
positions in vdB 106, vdB 133, and vdB 139 and at different nebular positions within vdB 
17 and vdB 59. We observe that /(5. 50-9. 75 /im)//(10. 25-14.0 /im) decreases by no more 
than 40% over Gq = 20 - 6 x 10^ 
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Our observational results are in strong agreement with Boulanger et al. (1998b), who 
concludes that the lEF flux ratios in ISO spectra of the Chamaeleon cloud, the p Oph 
cloud, NGC 7023, and M17 show no systematic dependence on Gq over Gq = 1 - 10^. 
Our results also strongly support Onaka (1999), who finds that the lEF flux ratios in ISO 
spectra of the Carina nebula are independent of Gq over Go = 100 - 10^. 

Much is yet to be gleaned from this ISO survey of reflection nebulae and especially the 
detection of lEFs toward vdB 133. In a future paper we will present a quantitative analysis 
of the lEF energy budget using laboratory absorption proflles of various carrier candidates 
(neutral PAHs, ionized PAHs, amorphous hydrocarbons, and coal tars). We will also 
present follow-up 7.75 /xm (LW6) ISO observations which cover an 8' by 8' mosaic around 
vdB 133, detailing the spatial distribution of the surrounding lEF emission. The relation 
of the lEF carriers to the gas and other components of dust surrounding vdB 133 will be 
explored by comparing the ISO images to ground-based infrared and molecular line (^^CO 
J =2-1) data that have also been obtained. 

We would like to thank A. Leger, Jean-Loup Puget and C Moutou for hosting us at 
the Institut d'Astrophysique Spatiale. We are grateful to W. Reach, D. Cesarsky and F. 
Boulanger for providing their expertise and invaluable help in the reduction of our ISO CAM 
data. We appreciate the generosity of D. Cesarsky and colleagues for providing us with 
their ISOCAM spectrum of NGC 7023, and F. Boulanger and colleagues for providing us 
with their ISOCAM spectrum of p Oph. 

We acknowledge NASA support of the ISO data analysis through NAG 5-3366 and 
JPL contract 961562, and also support from NATO Collaborative Research Grant 951347. 
This work was carried out in part at the Jet Propulsion Laboratory, California Institute of 
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Table 1. Stellar Data 



Nebula"^ 


Star 


Spectral Type 


TcB 


F(A<400nm) 
Ftotal 


F(A<240nm) 
Ftotal 


Aul„ 


{12Mm)'' 
(FIR) 


RcfsS 


Observed £ 


lources*^ 
















vdB17 


BD +30 549 


B8 V 


11,000 


0.55 


0.33 


0.09 




1 


vdB59 


HD 38563B 


B2 II-HI 


19,000 


0.87 


0.69 


0.14 




2 


vdB 101 


HD 146834 


G5HI 


5,000 


0.07 


0.003 


0.29 




1 


vdBlll 


HD 156697 


FOn 


7,300 


0.15 - 0.20^* 


0.02 - 0.03 


0.21 




1 


vdB 133 


HD 195593(A + B) 


F5 lab + B7 H 


6,800 + 12,000 


0.22'= 


0.09"= 


0.19 




1,3 


vdB 135 


BD +31 4152 


Ml Hie 


3,600 


0.01 


0.00005 


0.11 




1 


Comparison Sources 
















vdB 139 


HD 200775 


B3IVe 


17,000 


0.84 


0.65 


0.09 




1,4 


vdB 106 


HD 147933 + HD 147934 


B2 IV-V + B2 IV-V 


22,000 + 22,000 


0.92<= 


0.77^ 






5 



'^vdB 17 = NGC 1333, vdB 59 = NGC 2068, vdB 106 =pOph, and vdB 139 = NGC 7023. 

''Based on data from Tables 2 and 3 in ScUgrcn ct al. 1990. The adopted width of the 12 ^m IRAS bandpass is 1.75 X 10^'' Hz. 

'^ Sources observed in this study 

■^The range in Ft x<400nin) / ^total a-nd F/ x<240nin) / Ftotal stems from the uncertainty in the luminosity class of vdB 111. 

° F( XKiOOnm) / Ftotal for vdB 133 is equal to the sum of i^(A<400nm) for stars A and B divided by the sum of F^^tal for stars A and B. 
F(\<i2-l(inra) / Ftotal for vdB 133 is equal to the sum of i^(A<240nm) for stars A and B divided by the sum of F^^tal for stars A and B. 

Ft x<4W)nin) / Ftotal for p Oph is equal to the sum of i^(A<400nm) for stars HD 147933 and HD 147934 divided by the sum of Ff^tal for stars 
HD 147933 and HD 147934. F(A<240nm) /^totai for p Oph is equal to the sum of F(A<240nm) for stars HD 147933 and HD 147934 divided by 
the sum of Ftotal for stars HD 147933 and HD 147934. 

^Spectral classification references: (1) Racine 1968; (2) Strom et al. 1975; (3) Uchida, Sellgren, & Werner 1998 (Paper I); 
(4) Witt & Cottrell 1980; (5) de Geus, de Zeeuw, & Lub 1989. 
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Table 2. Observation Sequences 



Sequence 



Wavelengths (/xm)" 



Observed Sources (vdB) 



la 
lb 
2 
3 

4 



5.8, 6.2, 7.0, 9.5 

4.5^ 6.2, 7.0, 9.5 

7.0, 7.6, 8.4 

8.4, 8.62, 9.5 

10.5, 11.22, 12.0, 15.0^ 



101, 133, 135 

111 

101, 111, 133, 135 

101, 111, 133 
101, 111, 133, 135 



^Unless otherwise noted all wavelengths refer to circular variable filter settings, which 
have a spectral resolution A/AA = 40. 

^Narrowband filter LWl (AA = 1 fim). 
■^Narrowband filter LW9 (AA = 2 /im). 



Table 3. Field Positions 



Nebula 


Stellar Positioii 


I 


Source 


offset from star 


Sky offset from star 




RA (2000) 


Dec (2000) 


ARA( 


") 


ADec ( " ) 


ARA ( " ) 


ADec ( " ) 


Position Switched Observations 














vdB 101 


16 19 07.5 


-20 13 03.9 










+562 


-112 


vdB 111 


17 18 52.7 


06 05 07.1 










+498 


+160 


vdB 133 


20 30 59.1 


36 56 09.2 










-1187 


+266 


vdB 135 


20 36 45.5 


32 27 17.3 










-1218 


-97 


Full CVF Scan Observations 














VdB 17 (NGC 1333) 


03 29 19.8 


31 24 56.0 


+12 




-26 


+493 


+12 


vdB 59 (NGC 2068) 


05 46 44.6 


00 05 22.6 


+92 




+12 


+1026 


+ 13 


vdB 101 


16 19 07.5 


-20 13 03.9 


+39 




+20 


+470 


-24 


vdB 133 


20 30 59.1 


36 56 09.2 


+134 




+67 


+704 


+ 102 


vdB 135 


20 36 45.5 


32 27 17.3 


-96 




+201 


-771 


-86 
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Table 4. The Integration Regions 



Nebula Region Center" 




Region Size 


PA*" 


Comments'^ 


ARA(") 


AD 


ec(") 


AX(") 


AY(") 


(°) 




Position Switched Observations 






150 


150 


9.4 




vdB 101 







ER, XS 


vdBlll 







150 


150 


85.9 


ER, XS 


vdB 133 +72 




+46 


66 


30 


20.7 


LE, XS 


vdB 135 







150 


150 


5.6 


ER, XS 


Full CVF Scan Observations 






60 


60 


77.7 




vdB 17 (NGC 1333) +14 




-22 


LE 


vdB 59 (NGC 2068) +90 




+28 


72 


78 


5.8 


LE 


vdB 101 +34 




+227 


192 


192 


10.1 


ER 


vdB 133 +85 




+43 


60 


48 


83.8 


LE 


vdB 135 -104 




+204 


192 


192 


65.5 


ER 



^ offset of the integration region center from the central star. 

^ The rotation of the Y-axis east of celestial north. 

'^ Describes the subregion used to produce the spectral averages: LE — containing localized 



lEF emission, ER — 
from the central star. 



the entire usable portion of the image, XS — excluding the emission 



Note. — Columns 2 and 3 give the location of the subregion center in ARA, ADec offsets 
from the illuminating star(s), columns 4 and 5 its AX and AY dimensions, respectively, and 
column 6 the position angle of its Y-axis with respect to celestial north. Column 7 provides 
the brief description of the region, whether it encompasses the entire usable area of the array 
(ER — entire region), or a localized patch of lEF emission (LE — localized emission). In 
those cases where the image contained the illuminating star (the position switched data), the 
star and the immediate region containing a spurious arc-like reflection feature from the star 
(typically pixels X,Y = [6—20, 11-24] of the array) were excluded (XS; excluded star) from 
the spatial integration. 
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Table 5. Nebular Illumination 



Nebula Nebula distance 9'^^^ Gq (91 nm < A < 240 nm)'' Kefs': 

(pc) (arcsec) 



vdB 17 (NGC 1333) 


500 


29 


175 


vdB 59 (NGC 2068) 


550 


93 


520 


vdB 133 


1500 


100 


75 


Observed Sources 









Observed Sources 

1 
2 
3, 4 

vdB 139 (NGC 7023) 440 43 1820 5 

vdB 106 (p Oph) 125 2400 40 6 

^Tlie angular distance between illuminating star and the observed nebular region. 

''Go, the UV (91 nm < A < 240 nm) radiation field at the nebula, due to the illuminating stars, in units of 
1.6 X 10^'^ ergs s^^ cm^^, the local interstellar radiation field as determined by Habing 1968 

" Distance references: (1) Strom, Grasdalen, & Strom 1974; (2) Racine 1968; (3) Humphreys 1978; (4) 
Garmany & Stencel 1992; (5) Whitcomb et al. 1981; (6) de Geus, de Zeeuw, & Lub 1989. 
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Table 6. Normalized Integrated Line Fluxes 



Observed Sources 



vdB 17" vdB 59 vdB 133 



Comparison Sources 



vdB 139 vdB 106 



Method 1 Fluxes'" 



lEFs 



Broad Bumps 



Continuum 



(6.2 ^m)/(11.3 ^m) 
(7.7 /im)/(11.3 fim) 
(8.6 /im)/(11.3 ^m) 
(12.7 Mm)/(11.3 nm) 
(6-9 /im)/(11.3 ^m) 
(11 - 13 /im)/(11.3 fim) 
(5.2 - 15.1 Atm)/(11.3 fim) 



1.4 


1.7 


1.0 


1.1 


1.2 


3.8 


4.4 


4.8 


3.4 


4.1 


2.0 


2.3 


2.9 


1.8 


2.0 


0.18 


0.30 


0.29 


0.21 


0.21 


3.0 


2.7 


2.6 


2.5 


2.7 


1.6 


1.2 


1.1 


1.4 


1.8 


6.1 


4.6 


7.1 


5.0 


6.5 



Method 2 Fluxes 



lEFs 



Continuum 



(6.2 fim)/{n.3 fim) 
(7.7 Atm)/(11.3 fim) 
(8.6 /im)/(11.3 ^lln) 
(12.7 ^J,m)/{11.3 urn) 
(5.2 - 15.1 Atm)/(11.3 Atm) 



1.2 


1.7 


1.3 


1.2 


1.4 


1.7 


2.0 


2.3 


2.0 


2.0 


0.28 


0.44 


0.29 


0.24 


0.36 


0.18 


0.18 


0.24 


0.30 


0.23 


14.0 


11.0 


15.0 


13.0 


15.0 



Method 3 Fluxes 



lEFs 



Continuum 



(6.2 /^m)/(11.3 ^lln) 


1.4 


1.8 


1.0 


1.2 


1.0 


(7.7 ^m)/(11.3 ^lm) 


2.3 


3.0 


2.4 


2.5 


2.1 


(8.6 Mm)/(11.3 ^lm) 


1.4 


1.7 


1.7 


1.3 


1.3 


(12.7 Mm)/(11.3 fim) 


0.87 


0.50 


0.50 


0.63 


0.60 


(5 - 15 Mm)/(11.3 ^m) 


3.6 


3.7 


2.3 


3.1 


2.6 



11.3 /im Line Fluxes'^ 



Method 1 (11.3 ^m) 
Method 2 (11.3 ^m) 
Method 3 (11.3 ^m) 



1.3C-3 


2.8c~3 


1.5c~4 


. . A 


3.3e-4 


1.3C-3 


3.2C-3 


1.5e-4 


. . A 


3.3e-4 


2.2c~3 


4.0C-3 


3.2e-4 


. . d 


6.9e-4 



Total Band Fluxes 



IRAS 12 ^lm<= 
ISO LW 2-f 



l.le-2 


2.0e-2 


1.4e-3 


. . .9 


3.0e-3 


7.6e-4 


1.7e-3 


9.6e-5 


. . .9 


1.9e-4 



^'vdB 17 = NGC 1333, vdB 59 = NGC 2068, vdB 139 = NGC 7023, and vdB 106 = p Oph. 

''Integrated line flux from listed spectral feature (ergs s^-'^ cm^^ sr^-*^) divided by the integrated flux within the 11.3 
/im feature. 

•^Total integrated 11.3 ^m line flux in units of ergs s^^ cm^^ sr^^. 

''No integrated 11.3 ^m line flux of vdB 139 is presented since only a scaled version of the spectrum v^as available to 
us. 

''Total integrated source flux (ergs s^' cm^^ sr^') expected in 12 /^m IRAS band, found by the convolution of the 
source CVF spectrum vifith the 12 /im IRAS bandpass. 

'Total integrated source flux (ergs s^' cm~^ sr^') expected in ISO LW2 band, found by the convolution of the 
source CVF spectrum with the ISO LW2 bandpass. 

^No integrated 12 ^m line flux of vdB 139 is presented since only a scaled version of the spectrum was available to 
us. 
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Table 7. Normalized Bandpass Fluxes 



IRAS 12 ^m 



Filter 



ISO LW 2 



Method 1 Fluxes 

IEFs/(TotaP) 
Broad Bumps/Total 
Continuum/Total 

Method 2 Fluxes 

lEFs/Total 
Continuum/Total 

Method 3 Fluxes 

lEFs/Total 
Continuum/Total 



0.35 - 0.51" 
0.10-0.18 
0.33 - 0.45 



0.18 - 0.26 
0.74 - 0.82 



0.58 - 0.67 
0.29 - 0.47 



0.52 - 0.59 
0.22 - 0.26 
0.18 - 0.22 



0.26 - 0.35 
0.65 - 0.73 



0.64 - 0.71 
0.15 - 0.22 



'^Integrated flux of the spectral feature (erg cm^ s^ sr^ ) divided by the 
integrated flux within the broadband filter (IRAS 12 /im or ISO LW2) found by 
the convolution of the source CVF spectrum with the bandpass function. 

*^The range of values presented represent the minimum and maximum from 
among the five sources, vdB17 (NGC1333), vdB59 (NGC2068), vdB 133, 
vdB139 (NGC7023), vdB 106 (p Oph). 
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Flux time sequence (Seq 3) toward vdB 133 



X) 

o 
o 




500 1000 

Frame Number 



1500 



2000 



Fig. 1. — Rectification of ISOCAM data for Sequence 3 for vdB 133, during a 2 cycle 
sky -^ source observational sequence done consecutively with three different circular variable 
filter (CVF) settings at wavelengths of 8.6, 8.4, and 9.6 /im. The horizontal axis is image 
frame number; each image was obtained with an integration time of 0.28 s. a: the ideal 
response of the detector, b: an actual observation sequence toward vdB 133 showing the 
deviation from the ideal due to the detector flux transient response, c: the data of panel b 
after application of the flux "rectification" procedure. 
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Wavelength {/J-m) 



Fig. 2. — Flux-rectified and sky-subtracted 5.1-15 /im spectra of vdB 133, in Ii, (MJy sr~^) 
vs. A(/ini). Wavelengtlis of the 6.2, 7.7, 8.6, 11.3, and 12.7 /im infrared emission features 
(lEFs) are marked {dotted vertical lines), (a) Spectrum from position-switched observations 
at chosen CVF wavelengths {filled circles) . Error bars show the difference between the cycle 1 
and cycle 2 source— sky rectified data at each wavelength, {b) Spectrum from complete CVF 
scan {solid line). The difference between the rectified increasing and decreasing wavelength 
scans of the source-sky spectra {dotted line) is plotted as an estimate of the uncertainty. 
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10 12 

Wavelength (/Otm) 



16 



Fig. 3. — Flux-rectified, sky-subtracted 5.1-15 /xm spectra {solid line) from complete CVF 



scans of (a) vdB 17 (NGC1333), and (b) vdB59 (NGC2068), in /^ (MJy sr" 



vs. 



A(/im). 



The wavelengths of the 6.2, 7.7, 8.6, 11.3, and 12.7 /xm lEFs are marked {dotted vertical 
lines). The difference between the rectified increasing and decreasing wavelength scans of 
the source-sky spectra {dotted line) is plotted as an estimate of the uncertainty. 
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(b) vdB 101 (CVF) 




(d) vdB 135 (CVF) 




6 8 10 12 14 16 

Wavelength (/xm) 



4 6 8 10 12 14 16 

Wavelength (/i,m) 



Fig. 4. — Flux-rectified and sky-subtracted 5.1 - 15 fira spectra of {a, b) vdB 101, (c, 
d) vdB 135, and (e) vdB 111, in 1,^ (MJy sr^-*^) vs. A(/im). Wavelengtlis of the 6.2, 7.7, 
8.6, 11.3, and 12.7 ^m. lEFs are marked {dotted vertical lines), {a, c, e) Spectra from 
position-switclied observations at chosen CVF wavelengths {filled circles). Error bars show 
the difference between the cycle 1 and cycle 2 source— sky rectified data at each wavelength. 
{b, d) Spectra from complete CVF scan {solid line). The difference between the rectified 
increasing and decreasing wavelength scans of the source-sky spectra {dotted line) is plotted 
as an estimate of the uncertainty. 
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Fig. 5. — Plots of the fraction of UV emission to total emission as a function of stellar 
temperature, calculated from {^connected dots) Kurucz (1979) atmospheric models with log{g) 
= 4.0 and from {solid line) a blackbody. The UV range is defined as (a) A < 400 nm or (6) 
A < 240 nm. 
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Fig. 6. — Superposed spectra of [dotted line) vdB 17 (NGC 1333), (thick solid line) vdB 59 
(NGC 2068), and {thin solid line) vdB 133. Spectra are I^ (MJy sr^^), divided by their 
intensity in the IRAS 12 /im broad-band filter (MJy sr^^), vs. A(/im). The fraction of total 
stellar flux emitted at UV wavelengths, for the illuminating stars of these reflection nebulae, 
is lowest for vdB 133 and highest for vdB 59. 
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Model fits to vdb 17 (NGC 1333) 
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Fig. 7. — The 5-15 /iin spectrum of vdB 17 (NGC 1333) superposed with feature and 
continuum components fitted by (a) Method 1, {b) Method 2, and (c) Method 3 (see text). 
The integration regions adopted for the various lEFs in Method 3 (c) are marked {horizontal 
bars). The ISO LW2 and IRAS 12 /im filter transmission functions are also plotted {d). Note 
spectra are plotted as /A(erg cm~^ s~^ sr~^) vs. A(/xm), which differs from other figures. 
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Fig. 8. — Plot displaying log[-F(8.6 fim.)/F{11.3 /um)], the ratio of integrated 8.6 /im and 11.3 
/im lEF fluxes (erg cm^^ s^^ sr^^), vs. log(Tcfj), the effective temperature of the illuminating 
star(s) of each reflection nebula. Values of log[F(8.6 yum)/F(11.3 /im)] determined by Method 
1 {filled squares), Method 2 {filled triangles), and Method 3 {open circles) are compared (see 
text). Values of log[-F(8.6 fim.)/F{11.3 /im)] for each Method are connected {solid lines) for 
clarity. Values of log[F(8.6 /im)/F(11.3 /xm)] are plotted at two values of Tcfj, corresponding 
to the hotter and cooler illuminating star of vdB 133 (points connected by dotted lines). 
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vdb 17 (Radial Bins 



vdb 59 (Radial Bins) 




7 9 11 13 15 

Wavelength (/xm) 



Fig. 9. — Spectra {solid lines) of vdB 17 (left) and vdB59 (right), averaged over radial 
bins of projected distance from the illuminating star. Radial bins are characterized by the 
UV field, Go, incident on each nebular annulus (see text). Go is in units of 1.6 x 10~^ ergs 
s^^ cm~^, the local interstellar radiation field as determined by Habing (1968). Bins with 
values of Go = 20, 60, 200, 600, 1500, and 4000 (top to bottom) are shown for vdB 17. Bins 
with values of Go = 200, 600, 1700, 5000, 2.0 x 10^ and 6.0 x 10^ (top to bottom) are 
shown for vdB 59. Spectra are I^, (MJy sr^^), divided by the peak value of I^ at 7, 8.6, 11.3, 
and 12.7 /xm lEFs are marked [dotted vertical lines). The difference between the rectified 
increasing and decreasing wavelength scans of the source-sky spectra [dotted line) is plotted 
as an estimate of the uncertainty. 
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Fig. 10. — Full-width at half-maximuin (FWHM) of the 6.2 /xm, 7.7 /iin, and 8.6 /iin lEFs 
vs. incident UV intensity, log(Go); for different nebular positions within vdB17. FWHM 
were derived by Lorentzian fits to the spectra shown in Fig. ^ (see text). Error bars for each 
radial bin or annulus are derived from the difference between spectra of the top half of the 
annulus and the bottom half of the annulus. 
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Fig. 11. — Ratio of the integrated spectral intensity at 5.50 - 9.75 /xm to the integrated 
spectral intensity at 10.25 - 14.0 /im, 1(5.50-9.75 /im)//(10. 25-14.0 /im), vs. the incident 
UV intensity, log(G'o)- Ratios at different nebular positions within vdB 17 (NGC 1333; filled 
diamonds) and vdB59 (NGC 2068; filled circles) are derived from spectra in Fig. ^ Ratios 
for the spatially averaged spectra of vdB 133 {open circle; Fig. ^6 and Paper I), vdB 106 = 
p Oph [open square; Boulanger et al. 1996), and vdB 139 = NGC 7023 {cross; D. Cesarsky 
et al. 1996) are also plotted. Error bars for vdB 133 are derived from the difference between 
the rectified increasing and decreasing wavelength scans of the source-sky spectra. Error 
bars for each radial bin or annulus in vdB 17 are derived from the difference between spectra 
of the top half of the annulus and the bottom half of the annulus. Error bars for each radial 
bin or annulus in vdB 59 are derived from the difference between spectra of the right half of 
the annulus and the left half of the annulus. 



